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In 2008, researchers at HP Labs published a paper in Nature reporting the realisation of 
a new basic circuit element that completes the missing link between charge and flux- 
linkage, which was postulated by Leon Chua in 1971. The HP memristor is based on a 
nanometer scale Ti02 thin-film, containing a doped region and an undoped region. Fur- 
ther to proposed applications of memristors in artificial biological systems and nonvolatile 
RAM (NVRAM), they also enable reconfigurable nanoelectronics. Moreover, memristors 
provide new paradigms in application specific integrated circuits (ASICs) and field pro- 
grammable gate arrays (FPGAs). A significant reduction in area with an unprecedented 
memory capacity and device density are the potential advantages of memristors for In- 
tegrated Circuits (ICs). This work reviews the memristor and provides mathematical and 
SPICE models for memristors. Insight into the memristor device is given via recalling the 
quasi-static expansion of Maxwell's equations. We also review Chua's arguments based on 
electromagnetic theory. 
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1. Introduction 



Based on the International Technology Roadmap for Semiconductors (ITRS) report [ITRS 



200711 . it is predicted that by 2019, 16 nm half-pitch Dynamic Random Access Mem- 
ory (DRAM) cells will provide a capacity around 46 GB/cm 2 , assuming 100% area effi- 
ciency. Interestingly, memristors promise extremely high capacity more t han 110 GB /cm 2 
and 460 GB/cm 2 for 10 nm and 5 nm half -pitch devices, respectively IWilliamsl 120081 



Lewis and Lee , l2009tl . In contrast to DRAM memory, memristors provide nonvolatile op- 



eration as is the case for flash memories. Hence, such devices can continue the legacy 
of Moore's law for another decade. Furthermore, inclusion of molecular electronics and 
computing as an alternative to CMOS tech nologies in t he recent IT RS report emphasises 



the significant challenges of device scaling [Jones, 2009]. Moreover. lSwaroop et alJ H 1998 1 
demonstrated that the complexity of a synapse, in an analog VLSI neural network imple- 
mentation, is minimised by using a device called the Programmable Metallization Cell 
(PMC). This is an ionic programmable resistive device and a memristor can be employed 
to play the same role. 
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Research on memristor applications in various areas of circuit design, alternative mate- 
rials and spintronic memristors, and especially memristor device/circuit modeling have ap- 
peared i n the recent literature, (i) SPICE macro-modeling using linear and non-linear drift 



mode ls lichen and WangLl2009llZhang et all 120091. iBiolek et al.Ll2009bllBenderli and Wev 
2009 , Kavehei et al. , 2009 , Wang et all l2009all~ ii') Ap plication of memristors in differ- 
ent circuit configurations and their dy namic behaviour iLi et al. . I2009L ISun et al.l ■ l2009bl 



Wang et all l2009bi ISun et al 



200 9a|], (iii) Applic ation of memristor based dynamic sys- 



Lin and Wangl ll2Q09tl . (iv) Fine resolution programmable re 



Shin et al. ll2009 |, ( v) Me mrist or-based op - amp ci rcuit and 



tems to image encryption in 
sistor using a memristor in 
filter characteristics of memristors by lYu et al] II2009H and iWang et al] I2009cll . respec 
lively , (vi) Memristor receiver (MRX) stru cture for ultra-wide band (UWB) wireless sys 
tems I Itoh and Chual 12008 , Witrisal , 2009ll . (vii) Memristive system I/O nonlinearity can 



Delgadoli2009ll. (viii) The n umber of required memristors to compute a / : 
Lehtonen and Laiho l j 2009 1, and its digital log ic implementation using a memristor- 
ilRi 



cellation in 

function bv lLentonen ana Lainoillzuu ^|| . ana its dig 

based crossbar architecture in iRaia and Mo urad I20 09J1, (ix) Different phy s ical me cha 
nisms to store information in memristors IDriscoll et al. , I2009L IWang et all l2009dll . (x) 



Interesting fabricated nonvolatile memory using a. flexible memris tor, which is an inex 



pensi ve and low-power device solution is also recently reported BGergel-Hackett et al 



2009h . (xi) Using the memductance conce pt to develop an equivalent circuit diagram of 



a transmission-line has been carried out by Mallegoll 1 2003 . appendix 2] 



There are still many problems associated with memristor device level. For instance, it 
is not clear that which leakage mechanisms are associated with the device. The flexible 
memristor is able to retain its nonvolatile char acteristic for up to 14 days or, equivalently, 
up to 4000 flexes I Gergel-Hackett et al. . 2009tl. Thus, the nonvolatility feature eventually 



vanishes after a short period. I Strukov and Williams l2009ll also investigated this particular 



feature as a ratio of volatility to switching time. 

One disadvantage of using memristors is switching speed. The volati lity-to-switching 



speed ratio for memristor cells, in the HP cross-bar structure is aro und 10 II Strukov and Williams , 



2OO9I1 . which is much lower than the ratio for DRAM cells, 10 6 ilTRSll2007ll . therefore 



the switching speed of memristors is far behind DRAM. However, unlike DRAM, RRAM 



is non -volatile. In terms of yield, DRAM and RRAM are almost equal ILewis and Lee , 
Generally, endurance becomes very important once we note that the DRAM cells 
must be refreshed at l east e very 16 ms, which means at least 10 10 write cycles in their life 



time BLewis and Lee , l2009l l. Unfortunately, memristors are far behind DRAM in terms of 



the endurance view point, but the H P team is confide nt that there is no functional limitation 



against improvement of memristor I Williams. 1200811 . Finally, another advantage of RRAM 



is readability. Readability refers to the ability of a memory cell to report its state. This 
noise immunity in DRAM is weak because each DRAM cell stores a very small amount 
of charge particularly at nanometer dimensions, while in RRAM c ells, e.g. memristors, i t 
depends on the difference between the on and off state resi stances ]Lewis and Lee , 2009ll . 



This difference was reported up to one order of magnitude llWimamsLl200il 



e.g. [Buot and Raiago 


pal. 1994 Becketal. 


, 2000 


, Kriegeretall 2001 


Waser and Masakazul 


2007T Ignatiev et all 


2008[ 


Tulinaetal.. 2008ll. t 



tists were the only group that found the link between their work and the missing memris- 
tor postulated by Chua. Moreover, it should be noted that physically realised memristors 
must meet the math ematical requ ire ments of a memr i stor d evice or memristive systems 
that are discussed in Kang 1 1975 1 and Chua and Kang 1 1976ll . For instance, the hysteresis 
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loop m ust have a double- valued bow-tie trajectory. However, for example in lKrieger et al 



|2001] , the hysteresis loop shows more than two values for some applied voltage val- 
ues. Beck et al. I l2000ll demonstrates one of the perfect examples of memristor devices. iKim et al 



I 2009all recently introduced a multilevel one-time programmable (OTP) oxide diode for 
crossbar memories. They focused on a one-time programmable structure that basically 
utilises one diode and one resistor, 1D-1R, since obtaining a stable device for handling 
multiple programming and erasing processes is much more difficult than a one-time pro- 
grammable device. In terms of functionality, OTP devices are very similar to memristive 
elements, but in terms of flexibility, memristors are able to handle multiple programming 
and erasing processes. 

This paper focuses on the memristor device and reviews its device le vel properties . 
Although the memristor as a device is new, it was conceptually postulated bv lChual 1197111 . 
Chua predicted that a memristor could be realised as a purely dissipative device as a fourth 
fundamental circuit element. Thirty seven years later, Stan Williams and his group in the 
Information and Qua ntum Systems (IQS) Lab at HP realised the memristor in device form 
I Strukovet all l2008h . 

In Section |2j we review the memristor and its characteristics as a nano-switch, which 



was realised by Hewlett-Packard (HP) IStrukov et all 1200811 . and we review its properties 



based on the early mathematical models. We introduce a new model using a parameter we 
call the resistor modulation index (RMI). Due to the significance of ionic drift that plays 
the most important role in the memristive effect, this section is divided into two: ([aj» a lin- 
ear, and (0 a nonlinear drift model. Section [3]presents a preliminary SPICE macro-model 
of the memristor and different types of circuit elements in combination with a proposed 
memristor macro-model. Section [4] describes an interpretation of the memristor based on 
electromagnetic theory by recalling Maxwell's equations. Finally, we summarise this re- 
view in Section[5] 



2. Memristor Device Properties 



Traditionally there are only three fundamental passive circuit elements: capacitors, resis- 
tors, and inductors, discovered in 1745, 1827, and 1831, respectively. However, one can 
set up five different mathematical relations between the four fundamental circuit variables: 
electric current ;, voltage v, electric charge q, and magnetic flux (p. For linear elements, 
/(v,i) = 0, f(v,q) = where = Cv), and f(i, f) = where v = f (<p = Li), 

indicate linear resistors, capacitors, and inductors, respectively. 

In 1971, Leon Chua, proposed that there should be a fourth fundamental passive circuit 
element to set up a mathematical relationship between q an d <p, which he named the mem- 
ristor (a portmanteau of memory and resistor) IIChuaL 1 1 97 111 . Chua predicted that a class of 
memristors might be realisable in the form of a pure solid-state device without an internal 
power supply. 

In 200 8, Williams et al., at Hewlett Packard, announced the first fabricated memris- 
tor device llStrukov et all 1200811 . However, a resistor with memory is not a new thing. If 



taking the example of nonvolatile memory, it dates bac k to 1960 when Ber nard Widrow 
introduced a new circuit element named the memistor IWidrow et al. . 1 196(1 The reason 



for choosing the name of memistor is exactly the same as the memristor, a resistor with 
memory. The memistor has three terminals and its resistance is controlled by the time inte- 
gral of a control current signal. This means that the resistance of the memistor is controlled 
by charge. Widrow devised the memistor as an electrolytic memory element to form a 
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basic structure for a neural circuit architecture called AD ALINE (ADAptive Linear NEu- 
ron), which was introduce d by him and his postgraduate student, Marcian Edward "Ted" 
Hoff I Widrow et al. . 1960ll . However, the memistor is not exactly what researchers were 
seeking at the nanoscale. It is just a charge-controlled three-terminal (transistor) device. 
In addition, a two-terminal nano-device can be fabricat ed without nanoscale align ment, 
which is an advantage over three-terminal nano-devices I Lehtonen and Laihol 2009 1 . Fur- 
thermore, the electrochemical memistors could not meet the requirement for the emerging 
trend of solid-state int egrated circuitry . 

Thirty years later, iThakoor et al. introduced a solid-state thin-film tungsten- 

oxide-based, nonvolatile memory. The concept is almost similar to the HP memristors. 
Their solid-state memistor is electrically reprogrammable, it has variable resistance, and it 
is an analog synaptic memory connection that can be used in electrical neural networks. 
They claimed that the resistance of the device could be stabilised at any value, between 
100 k£2 and 1 G£2. This solid-state memistor has a thick, 60-80 nm, layer of silicon diox- 
ide that achieves nonvolatility over a period of several months. This thick electron blocking 
layer, however, causes a considerable reduction in the applied electric field. As a conse- 
quence, they reported very high programming voltages (± 25 to 30 ) and very long (min- 
utes to hours) switching times. 

In the 1960s, the very first rep ort on the hysteresis behaviour of current- voltage curve 
was published bv lSimmonsI 111 96311 . which is known as the Simmons tunneling theory. The 

Simmons theory generally characterises the tunneling current in a Meta l-Insulator-Metal 

(MIM) system. A variable resistance with hysteresis was also published bv lSimmons and Verderber 
1 196711 . They introduced a thin-film (20 to 300 nm) silicon dioxide doped with gold ions 
sandwiched between two 200 nm metal contacts. Thus overall the device is a MIM system, 
using aluminum metal contacts. Inte restingly, their sy stem demonstrates a sinh function 
behaviour as also recently reported in lYang et al.1 H2008I1 for the HP memristor. It is also re- 
ported that the switching from high- to low impedance takes about 100 ns. As their device 
is modelled as an energy storage element it cannot be a memristor, becau se memristor s 
remember the total charge passing through the port and do not store charge lOsteri 1 1 974 1. 

The sinh resistance b ehaviour of memristors can b e utilised to compensate the linear- 
ity of analog circuits. In IVarghese and Gan dhi Il20q9|j a memristor based amplifier was 
proposed utilising th e behaviour of a sinh resistor jTavakoli and Sarpeshkai , 2005 1 as a 

memristor element. In Sh in et al. |2009| such a structure was also introduce d as an example 

withou t mentioning the sinh resistance behaviour of memristors. The idea in lTavakoli and Sarpeshkai 
Il2005ll is to characterise a sinh function type circuit that can be used to linearise a tanh 
function type circuit behaviour, e.g. CMOS differential amplifier. A theoretical analysis 
shows that a sinh function significantly reduces the third harmonic coefficient and as a 
consequence reduces nonlinearity of circuit. 

Chua mathematically predicted that there is a solid-sta te device, which defines the miss- 
ing relationship between four basic variables | Chual 1971 ] . Recall that a resistor establishes 
a relation between voltage and current, a capacitor establishes a charge-voltage relation, 
and an inductor realises a current-flux relationship, as illustrated in Fig. Q] Notice that we 
are specifically discussing nonlinear circuit elements here. 

Consequently, <p = /m(<?) or q = gM(<P) defines a charge-controlled (flux-controlled) 



memristor. Then, 



dt 



f£/k(gj dq 
dq 



dq _ dg M {<p) d<p 
d<p dt 



, which implies v(f) 



K*) = dS d!p l ' ) v ( t )- Note mat ' M (<7) = ^l!r for a charge-controlled memristor and W((p) = 



^2q> f° r a flux-controlled memristor, where M(q) is the incremental memristance and 



Article submitted to Royal Society 



The Fourth Element 



5 



Figure 1. The four fundamental two-terminal circuit elements. There are six possible relationships 
between the four fundamental circuit variables, current (£), voltage (v), charge (q), and magnetic flux 
(<p). Out of these, five relationships are known comprising, i = t? and v = and three circuit 
elements: resistor (R), capacitor (C), and inductor (L). Resistance is the rate of voltage change with 
current, inductance is the rate of flux change with current, capacitance is the rate of charge change 
with voltage. The only missing relationship is the rate o f flux change with vo ltage (or vice- versa) that 
leads to the fourth element, a memristor. Adapted from lStrukov et al] J2008h . 



W(<p) is the increment al memducta nce where M(q) is in the units of Ohms and the units 
of W(<p) is in Siemens llChuall 197 1 1 . 

Memristance, M(q), is the slope of the (p-q curve, which for a simple case, is a piece- 
wise (p-q curve with two different slopes. Thus, there are two different values for M(q), 
which is exactly what is needed in binary l ogic. F or detailed information regarding typical 
(p-q curves, the reader is referred to lChual 11971 1. 

It is also obvious that if M{q) > 0, then instantaneous power dissipated by a memristor, 
p(i) = M(q)(i{t)) 2 , is always positive so the memristor is a passive device. Thus the (p-q 
curve is a monotonically-i ncreasing function. T his feature is exactly what is observed in 



the HP memristor device flStrukov et all 1200811 . Some other properties of t he memristor 



such as zero-crossing betw een current and voltage signals can be found in IChual 11971 ] 



sing betw 

and lChua and Kana jl976ll . T he most important featu re of a memristor is its pinched hys- 
teresis loop v-i characteristic I Chua and Kang , 1976ll . A very simple consequence of this 



property and M(q) > is that such a device is purely dissipative like a resistor. 

Another important property of a memristor is its excitation frequency. It has been 



show n that the pinched hysteresis loop is shrunk by increasing the excitation frequency IChua and K ang 



197611 . I n fact, when the freque ncy increases toward infinity, the memristor acts as a linear 
resistor IChua and Kangll 197611 . 



Interestingly enough, an attractive property of the HP memristor IStrukov et all 1200811 



which is exclusively based on its fabrication process, c an be deduced from the HP mem- 
ristor simple mathematical model llStrukov et all 120081 and is given by, 



M{q)=R OFF [ 1- 



(i ->.)), 



(2.1) 



where j3 has the dimensions of magnetic flux (p(t). Here, j3 = ^ in units of sV = Wb 



where /Zd is the average drift mobility in unit of cm /sV and D is the film (titanium 
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dioxide, TiC>2) thickness. Note that /?off an d ^on are simply the 'on' and 'off state re- 
sistances as indicated in Fig. [2] Also q(t) defines the total charge passing through the 
mem ristor device in a time window, t - to. Notice tha t the memris tor has an internal 
state BChua and Kangl 1197611 . Furthermore, as stated in lOsterl 1U974I1 . q{t) = J t(j i(T)dT, 
as the state variable in a charge-controlled memristor gives the charge passing through 
the device and does not behave as storage charge a s in a capacitor as incorrectly reported 
in some works, e.g. lSimmons and Verderberl II 196711 . This concept is very important from 
two points of view. First of all, a memristor is not an energy-storage element. Second, this 
shows that the memristor is not merely a nonlinear resistor, but is a nonlinear resistor with 



charge as a state variable lOsten, 1197411 



Five years after Chua's paper on the memristor IIChualll971ll . he and his graduate stu- 
dent, Kang, published a paper defining a much broader class of systems, named memristive 

systems. From the memristive systems viewpoint a generalized definition of a m emristor is 

v(t ) = R(w)i(t), where w defines the internal state ofthe system and ^ = f(w,i) flChua and Kang . 
Based on this de finition a mem r istor i s a special case of a memristive system. 



The HP memristor IStrukov et all 1200811 can be defined in terms of memristive sys- 



tems. It exploits a very thin-film TiC>2 sandwiched between two platinum (Pt) contacts and 
one side of the TiC>2 is doped with oxygen vacancies, which are positively charged ions. 
Therefore, there is a TiC>2 junction where one side is doped and the other is undoped. Such 
a doping process results in two different resistances: one is a high resistance (undoped) and 
the other is a low resistance (doped). Hence, HP intentionally established a device that is 
illustrated in Fig. [2] The internal state variable, w, is also the variable length of the doped 
region. Roughly speaking, when iv->0we have nonconductive channel and when w —> D 
w e have conductive channel. The HP memristor switching mechanism is further discussed 
in 



Yang et all 120081. 



W 



D 



Pt 


Doped 


Undoped 


Pt 



Figure 2. Schematic of HP memristor, where D is the device channel length and w is the length 
of doped region. The size of doped region is a function of the applied charge and is responsible 
for memristive effect as it changes the effective resistance of the dev ice. Usually w is shown by its 
normalised counterpart, x = w/D. Adapted from lStrukov et all l2008ll . 



(a) Linear drift model 



T he memristor's state equation is at the heart of the HP memristive system mecha- 
nism I Wang , 2008 , Joglekar and Wolf , 2009ll . Let us assume a uniform electric field across 
the device, t hus there is a linear relatio nship between drift-diffusion velocity and the net 
electric field IBlanc and Staeblen.1197111 . Therefore the state equation is, 



1 dw(t) _ Ron 
D dt ~ B 



(2.2) 



Article submitted to Royal Society 



The Fourth Element 



7 



Integrating Eq. 12.21 gives = + ^jf-q{t), where w(fo) is the initial length of 
w. Hence, the speed of drift under a uniform electric field across the device is given by 
\>d = ■ In a uniform field we have D = Vd x f. In this case 2d = ixt also defines the 
amount of required charge to move the boundary from w(?o), where w — > 0, to distance 
w(ti)), where w — > D. Therefore, Qd = ^j^-, so 

w(t) _ w(t ) , <?(/) 



If jc(f) = ^ then 



D D Q D 



x(t)=x(t ) + ^-, (2.4) 
2d 



where ^ describes the amount of charge that is passed through the channel over the 
required c harge for a conductiv e channel. 



Using IStrukov et al.l 1200811 we have, 



D 



v{t) = [Ron ^T+RoffI 1 - I K0 ■ (2-5) 



By inserting a: (f) = Eq. l2.5l can be rewritten as 

v(t) = \ R Kx(t)+R 0¥P (l -x(t)X\i(t) . (2.6) 



Now assume that q(to) = then w(t ) = w(/o) 7^ 0, and from Eq. 12.61 

M =/?oN^(fo) + r^l-x(f )^ , (2.7) 

where r = and Mo is the memristance value at to. Consequently, the following equation 
gives the memristance at time t , 

M(q)=Mo-AR\^J , (2.8) 



where AR = Roff — Ron- When ^off ^ Ron. Mq ~ ^off an d Eq. l2.1l can be derived 
from Eq.EJ 

Substituting Eq.^into v(f ) = M{q)i(t), when i(t ) = ^f-, we have, 
Recalling that M(q) = d y q ^ , the solution is 



/ s QdM / / 2AR , s \ 

4 f = ^— Si 1 ± J I -m( t ) . (2.10) 

iy> AR \ V QdMq J 
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Using AR w Mq /?off> Eq. I2.10l becomes. 

rf'-^'-^-db*®)- (2U) 

Consequently, using Eq. 12 .41 if Qo = ^j^ qn , so the internal state of the memristor is 




x(f) = l-Ul-^ 9 (f) ] . (2.12) 

The current-voltage relationship in this case is, 

v(t) 

iff) = 7 ^ 

^OFF 

In Eq. 12.131 the inverse square relationship between memristance and Ti02 thickness, 
D, shows that for smaller values of D, the memristance shows improved characteristics, 
and because of this reason the memristor imposes a small value of D. 

In Eqs. 12. 10112.131 the only term that significantly increases the role of <p(f) is lower 
Qo. This shows that at the micrometer scale Rq ^,q d = = j§i is negligible and the 
relationship between current and voltage reduces to a resistor equation. 

Substituting /3 = ^ that has the same units as magnetic flux into Eq. 12.131 and consid- 
<p{t) _ fi D <p{t 



ering c(t) = -^p- = ^ 2 as a normalised variable, we obtain 



v(t) 

i{t) = ^ r , (2.14) 



^OFF 



where y 1 — jc(t) is what we call the resistance modulation index, (RMI) iKavehei et al 
2009h . 



Due to the extremely uncertain nature of nanotechnologies, a variability-aware mod- 
eling approach should be always c onsidered. Two well-kn ow solutions to analyse a mem- 
ristive system were investigated in lChen and Wang l2009ll . 1) Monte-Carlo simulation for 



evaluating (almost) complete statistical behaviour of device, and 2) Corner analysis. Con- 
sidering the trade-off between time-complexity and accuracy between these two approach 
as shows the importance of using a simple and reasonably accurate model, because finding 
the real corners is highly dependent on the model accuracy. The resistance modulation in- 
dex, RMI, could be one of the device parameters in the model extraction phase, so it would 
hel p to provide a simp l e mod el with fewer parameters. 



ip i 
J< 



Jog lekar and W olf [2009] clarified the behaviour of two memristors in series. As shown 
in Fig. [3j they labeled the polarity of a memristor by 77 = ±1, where 77 = +1 signifies 
that w{t) increases with positive voltage or, in other words, the doped region in mem- 
ristor is expanding. If the doped region, w(t), is shrinking with positive voltage, then 
T) = — 1, In other words, reversing the voltage source terminals implies memristor po- 
larity switching. In Fig.[3](a), the doped regions are simultaneously shrunk so the overall 
memrist ive effect is retained. In F ig. [3](b), however, the overall memristive effect is sup- 



pressed ] Joglekar and WolfLl2009ll . 
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(a) o Doped 


Undoped ^H- 




(b) c ^| Doped 


Undoperi 



■^H Doped 


Undoped ^H. 




Undoped 


Doped 



Figure 3. Two memristors in series, (a) With the same polarity, both r\ = —1 or both ?l = +1. (b) 
With opposite polarities, rj = — 1 and TJ = +1. Where T] = +1 signifies that w(f) increases with 
positive voltage or, in other words, the doped region in memristor is expanding and 7] = — 1 indi- 
cates tha^_fiie_dorjed_reeionjs shrinking with applied positive voltage across the terminals. Adapted 
from lJoglekar and Wolf! l2009ll . 



(2.15) 



Using the memristor polarity effect and Eq. 12.21 we thus obtain 



1 dw(t) Ron .. . 



Then with similar approach we have 

m = - 



(2.16) 



floFF Jl-^c{t) 



There is also no phase shift between current and voltage signals, which implies that the 
hysteresis loop always crosses the origin as demonstrated in Fig. [4] For further investiga- 
tion, if a voltage, v(t) = vosin(of), is applied across the device, the magnetic flux would 
be (p{t) — — ^ cos(a)f). The inverse relation between flux and frequency shows that at very 
high frequencies there is only a linear resistor. 

Fig. IS de monstrates E q. l2~16l in MATLAB for five different frequencies, where a\) = 
l2008ll parameter values. D = 10 nm. jdn = 10" 10 cm 2 s" 1 V~ ! 
16 k£2, v'o = 1 V, and TJ = — 1. The resistance modulation index, 

RMI = 



using the lStrukov et al 
^on = 100 £2, Rqff 



1 — rjjc{t) simulation with the same parameter values is shown in Fig. [5] 
Using differe nt parameters ca uses a large difference in the hysteresis and memristor 
characteristics. In Witrisal ll2009tl . a new solution for ultra-wideband signals u s ing m em- 
ristor devices was introduced. Applying the parameter values given in lWitrisall j2009ll re- 
sults in a significant difference in the value of C0q. Substituting parameter values given 
in the paper gives (Oq w 4 GHz instead of On w 50 kHz, using to to rj 0.1 ms and the 
parameter values from IStrukov et al.l fl2008ll . The new parameter values are D = 3 nm, 



Hd = 3 x 10" 8 m 1 s" 1 V" 1 , Ron = 100 Q., R OFF = 10 k£2, and v = 0.2 V. Using these 
parameters shows that even though the highest and lowest memristance ratio in the last 
case, Fig.|U is around 2, here the ratio is approximately equal to 120. 



(b) Nonlinear drift model 

The electrical behaviour of the memristor is directly related to the shift in the boundary 
between doped and undoped regions, and the effectively variable total resistance of the 
device when an electrical field is applied. Basically, a few volts supply voltage across a 
very thin-film, e.g. 10 nm, causing a large electric field. For instance, it could be more than 
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Figure 4. The hysteresis of a memristor based on Eq. 12.161 This verifies the hysteresis shrinks at 
higher frequencies. It also shows that the effective resistance is changing, so there is a varying mem- 
ristance with a monotonically-increasing q-(p curve. 
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Figure 5. The hysteresis characteristics of the memristor. It shows that the memristance value is 
varying from a very low to a very high resistance. It is clear that these values depend on the parameter 
values, such as R on and R g, 



10 6 V/ cm, which results in a fast and significant reduction in energy barrier jBlanc and Staebleil 



2009] 



1971 ] . Therefore, it is reasonable to expect a high nonlinearity in ionic drift-diffusion I Waser et al 



One significant problem with the linear drift assumption is the boundaries. The linear 
drift model, Eq. 12.21 suffers from problematic boundary effects. Basically, at the nanoscale, 
a few vol ts causes a large ele ctric field that leads to a significant nonlinearity in ionic 
transport IStrukov et all 1200811 . A few attempts have been carried out so far to consider 
this nonlinearity in the state equation | Strukov et all 120081 IStrukov and Williams , I2009L 
Biolek et~all l2009bl Benderli and Wev , 2009 ] . All of them proposed a simple window 



function, F(|), which is multiplied by the right-hand side of Eq. 12.21 In general, £, could 
be a variable vector, e.g. £, = (w,i) where w and i are the memristor's state variable and 
current, respectively. 

In general, the window function can be multiplied by the right-hand side of the state 
variable equation, Eq. 12.21 
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(2.17) 



where x(t) = w(t) /D is the normalised form of the state variable. The window function 
makes a large difference between the model with linear and nonlinear drift assumptions at 
the boundaries. Fig. [6] shows such a condition considering a nonlinear drift assumption at 
the critical, or boundary, states. In other words, when a linear drift model is used, simula- 
tions should consider the boundaries and all constraints on initial current, initial voltage, 
maximum and minimum w, and etc. These constraints cause a large difference in output 
between linear and nonlinear drift assumptions. For example, it is impossible to achieve 
such a realistic curve, as in Fig. [6] using the linear drift modeling approach. 




c 

t 0.0 
u 




Figure 6. The hysteresis characteristics using the nonlinear drift assumption. This hysteresis shows 
a highly nonlinear relationship between current and voltage at the boundaries. To model this nonlin- 
earity, there is a need for an additional term on the right hand side of the memristor's state equation, 
called a window function. 



In Strukov et al 



j2008ll . the window function is a function of the state variable and it 
is defined as F(w) = w(l — w) /D . The boundary conditions in this case are F(Q) = and 
F(D) = i-j=p ss 0. It meets the essential boundary condition F(^ — > boundaries) = 0, except 
there is a slight difference when w — > D. The problem of this boundary assumption is when 
a memristor is driven to the terminal states, w — > and w — > D, then ^ — > 0, so no external 



field can change the memristor state IBiolek et al.L l2009bll . This is a fundamental problem 



of the window function. The second problem of the window function is it assumes that the 
memristor remembers the amount of charge that passed through the d evice. Basically, this 
is a direct result of the state equation, Eq. 12.21 IBiolek et all l2009bll . However, it seems 
that the device remembers the position of boundary between the two regions, and not the 

amou nt of charge. 

In lBenderli and Wevli2009ll. ano ther w indow function has been proposed that is slightly 
different from that in lStrukov et al.l 1200811 . This window function, F(w) = w(D — w)/D 2 , 
approaches zero when w — > and when w — > D then F(w) — > 0. Therefore, this window 
function meets both the boundary conditions. In fact, the second window function is imi- 
tates the first function when we consider x =w/D instead of w. In addition, there is another 
problem associated with these two window functions, namely, the modeling of approximate 
linear ionic drift when < w < D. Both of the window functions approximate nonlinear 
behaviour when the memristor is not in its terminal states, w = or w = D. This problem is 
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addressed in lJoglekar and Woln 11200911 where they propose an interesting window function 
to address the nonlinear and approximately linear ionic drift behaviour at the boundaries 
and when < w < D, respectively. Nonlinearity (or linearity) of their function can be con- 
trolled with a second parameter, which we call the control parameter, p. Their window 
function is F(x) = 1 — (2x — l) 2p , where x ~w/D and p is a positive integer. Fig. 0(a) 
demonstrates the function behaviour for different 2 < p < 10 values. This model considers 
a simple boundary condition, F(0) = F(l) = 0. As demonstrated, when p > 4, the state 
variable equation is an approximation of the linear drift assumption, F(0 <x < 1) ps 1. 

The most important problem associated with this model is revealed at the boundaries. 
Based on this model, when a memristor is at the terminal states, no external stimulus can 



change its state. iBiolek et al.l Il2009bll tackles this problem with a new window function 
that depends on x, p, and memristor current, i. Basically, x and p are playing the same role 
in their model and the only new idea is using current as an extra parameter. The window 
function is, F(x) = 1 — (jc— sgn(— i)) 2p , where i is memristor current and sgn(/) = 1 when 
i > 0, and sgn(/) = when ; < 0. As a matter of fact, when the current is positive, the 
doped region length, w, is expanding. Fig. 0(b) illustrates the window function behaviour. 




0.4 0.6 

x 

(a) (b) 

Figure 7. Non-linear window functions, (a) F(x) = 1 — (2jc— 1) 2 '\ (b) F(x) = 1 — (x— sgn(— i)) 2p . 
There are around four window functions in the literature but these two are meet the boundary condi- 
tions. 



All window functions suffer from a serious problem. They are only dependent on the 
state variable, x, which implies that the memristor remembers the entire charge that is 
passing through it. Moreover, based on the general definition of the time-invariant memris- 
tor's state equation and current-voltage relat ion, x = f(x, i) or x = g(x, v), / and g must be 
continuous n-dimensional vector functions I Kand, 19751 chap. 2]. However, the last win- 
dow function, F(x) = 1 — (x — sgn(— i)) 2p , does n ot provide the continuity condition at the 
boundaries, x — > or x — > 1 . IBiolek et al.l ll2009bll did not use the window function in their 
recent publication llBiolek et alll2009all . 

In Strukov and Williams! 120091 the overall drift velocity is identified with one linear 
equation and one highly nonlinear equation, v = jiE, when£ << £0 an d V = IxEoexp(^), 
for E ~ Eq, where V is the average drift velocity, E is an applied electric field, jj, is the 
mobility, and Eq is the characteristic field for a particular mobile ion. The v alue of Eo is 
typically about 1 MV/cm at room temperature I Strukov and Williams . 2009h . This equa- 
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tion shows that a very high electric field is needed for exponential ion transport. They also 
showed that the high electric field is lower than the critical field for dielectric breakdown 
of the oxide. Reviewing the available window functions indicates that there is a need for an 
appropriate model that can define memristor states for strongly nonlinear behavior where, 
E~E . 

In addition to the weakness of nonlinear modeling in the original HP model, there 
are some other drawbacks that show the connection between physics and electronic be- 
haviour was not well established. Moreover, the currently available electronic models for 
memristors followed the exact pathway of the HP modeling, which is mostly due to the 



fact th at the underlying physical conduction mechanism is not fully clear yet I Ki m et al 



2009bll . One weakness is that the HP model does not deliver any insight about capaci- 



tive effects, which are naturally associated with memristors. These capacitive effects will 
later be exp lained in terms of a memcapacitive effect in a class of circuit elements with 
memory. In lKim et al. I ll2009bll the memristor behaviour was realised using infinite num- 



ber of crystalline magnetite (FesO/t) nanoparticles. The device behaviour combines both 
memristive (time-varying resistance) and memcapacitive (time-varying capacitance) ef- 
fects, which deliver a better model for the nonlinear properties. Their model description 
for current-voltage relationship is given as, 

v(t) 

i(t) = U == , (2.18) 



where R(x,t) and C(x,t) are the time-varying resistance and capacitance effects, respec- 
tively. The time-dependent capacitor is a function of the state variable, x(t) = w(t)/D and 
AC(t ), where AC(t ) is defined as the additional capacitance caused by changing the value 
of capacitance llKim et al. , 2009bll . C(x,t) = ^y^^, where Con is the capacitance at 



x = 0. The state variable is also given by, 



Kim e t al. [ 2009b] also investigated the impact of temperature variation on their Fe304 
nanoparticle memristor assemblies of D equal to 9, 12, and 15 nm. It was reported that the 
change in electrical resistivity (specific electrical resistance), p r , as an explicit function of 
temperature can be defined by, logp, = 1 / \/T, which means there is a significantly in- 
creasing behaviour as temperature decreases. As a consequence, for example, there is no 
hysteresis loop signature at room temperature, T = 295 K, (D = 12 nm) while at T = 210 K 
it shows a nice bow-tie trajectory. As they claimed, the first room temperatur e reversible 



switch ing behaviour was observed in their nanoparticle memristive system IKim et al 
l2009bh . 

It is worth noting that there are also two other elements with memory named the mem- 
capacitor (Memcapac itive, MC, Systems) and meminductor (Meminductive, MI, Systems) 



Di Ventra et al.1 12009H postulated that these two elements also could be someday realised in 
device form. The main difference between these three elements, the memristor, memcapac- 
itor, and meminductor is that, the memristor is not a lossless memory device and dissipates 
energy as heat. However, at least in theory, t he memcapacitor and m eminductor are lossless 
devices because they do not have resistance. Di Ventra et al. I l2009ll also investigated some 
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examples of using different nanoparticle-based thin-film materials. Memristors (Memre- 
sistive, MR, Systems) are identified by a hysteresis current-voltage characteristic, whereas 
MC and MI systems introduce Lissajous curves for charge-voltage and flux-current, re- 
spectively. Similar to memristors, there are two types of these elements, therefore, the three 
circuit elements with memory (mem-systems/devices) can be summarised as follows, 

Memristors (MR Systems): A memristor is an one-port element whose instantaneous 
electric charge and flux-linkage, denoted by q mr and <p nir , respectively, satisfy the 
relatio n F mr (q mr , (p mr ) = 0- It has been proven that these devices are passive ele- 



ments OStrukov et all 1200811 . As discussed, they cannot store energy, so v mr (f) = 
whenever i mv (t) = and there is a pinched hysteresis loop between current and 
voltage. Thus, charge-flux curve is a monotonically-increasing function. A mem- 
ristor acts as a linear resistor when frequency goes toward infinity and as a non- 
linear resistor at low frequencies. Due to the nonlinear resistance effect, c ^ff- = 

R (t) + imi^W should be utilised instead of ^ = i?(/)^f . There are two types 
of control process, 

I. nth order current-controlled MR systems 0, q m r = Jimr(^)dT, 

• V mr (f) =R(x,imr, t)imr(t) 

• % ' = yicmr(-^ ^mrjO 

II. nth order voltage-controlled MR systems, <p ml - = / v ml {x)d%, 

• imr(t) =R~ 1 (x,V mr ,t)v mI (t) 

• X = /vcmrO^ Vmrj^J ■ 

Memcapacitors (MC Systems): A memcapacitor is an one-port element whose instan- 
taneous flux-linkage and time-integral of electric charge, denoted by <p mc and <7 mc , 
respectively, satisfy the relation F mc ((p mc , <7 mc ) = 0. The total added/removed energy 
to/from a voltage-controlled MC system, U mc = /v mc (T)j' mc (T)dT, is equal to the 
linear summation of areas between q mc -v mc curve in the first and third quadrants with 
opposite signs. Due to the nonlinear capacitance effect, = ; mc (f ) = C(*)^?p + 
Vmc(t)~fr should be utilised instead of ^ = C(t)^, so U mc = f \> mc Cdv mc + 
J v^ c dC. In prin ciple a memcapacitor can be a passive, an active, and even a dissi- 
pativ^l element [ Di Ventra et al. I l2009ll . If v mc (/) = V mc ocos(27TOf) and capacitance 



is varying between two constant values, Con and Coff, then the memcapacitor is 
a passive element. It is also important to note that, assuming zero charged initial 
state for a passive memcapacitor, th e amount of removed energy cannot exceed the 
amount of previously added energy IDi Ventra et al. L l2QQ9tl . A memcapacitor acts as 



a linear capacitor when frequency tend to infinity and as a nonlinear capacitor at low 
frequencies. There are two types of control process, 

I. nth order voltage-controlled MC systems, <p mc = / v mc (T)dT, 

• <?mc(0 = C(x,Vmc,t)v mc (t) 

• X = /vcmc (-^j ^mc? 

t C urrent- (and volt age-) controlled is a b etter definition for memristors because they do not store any 
charge IQster, IT974ll . In |Pi Ven tra et al J |2009tl it is specified as current- (and voltage-) controlled instead of 
charge- (and flux-) controlled in lChua tl 97 ill . 

$ Adding energy to system. 
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II. nth order charge-controlled MC systems, <J n 

• v mc (f) = C~ l (x,q mc ,t)q mc (t) 

* X = /qcmc (-^ j ?mc s ■ 



Jq mc (?)dT = JJi mc (?)dT, 



Meminductors (ML Systems): A meminductor is a one-port element whose instantaneous 
electric charge and time-integral of flux-linkage, denoted by q m \ and p m i, respec- 
tively, satisfy the relation F m \(p m \,q m \) = 0. In the total stored energy equation in a 
ML system, U m \ = J v m \{x)i m \{x)dx, the nonlinear inductive effect, = v m \{t) = 
L^JM. _|_ i^jyjL. s hould be taken into account. Thus, U m \ = Ji m \Ldi m \ + Ji^dL. 
Similar to MC systems, in principle a ML system can be a passive, an active, and 
even a dissipative element and usi ng the same app r oach, under some assumptions 
they behave like passive elements IIDi Ventra et all 1200911 . There are two types of 
control process, 

I. nth order current-controlled ML systems, q m \ = Ji m i(x)dT, 

• <jV(0 = L(x,i mU t)i m i(t) 

• * = /icml(*,jml,0 

II. nth order flux-controlled ML systems, p m \ = f (p m \(t)dT = ff v m i(T)dT, 

• ; m i(f) =ir 1 (x,<p ml ,t)<p ml (t) 

• -* = /fcml(*,Pml,*) • 



Biole k et al. I ll2009all introduced a generic SPICE model for mem-devices. Their mem- 



ristor model was discussed before in th is section. Unfortunately, there are no simulation 
results available in lBiolek et al.l l2009all for MC and ML systems. 

It is worth notin g that there is no equivalent mechanical element for the memristor. 
In Chen et al. i2009f a new mechanical suspension component, named a J-damper, has 
been studied. This new mechanical component has been introduced and tested in Formula 



One R acing, delivering significant performance gains in handling and grip IChen et al 



2009]. In that paper, the authors attempt to show t hat the , J-dam per, which was invented 
and used by the McLaren team, is in fact an inerter llSmilhll2002ll . The inerter is a one-port 
mechanical device where the equal and opposite applied force at the terminals is propor- 
tional to the relative acceleration between them I SmithL 2002 1 . Despite the fact that the 
missing mechanical circuit element h as been chosen because of a "spy scandal" in the 2007 
Formula One race IFormula ll 1200711 . it may be possible that a "real" missing mechanical 
eq uivalent to a me mristor may s omeday be found, as its mechanical model is described 
bv lOsted lll973h and lOsterl Il974 l. 



3. SPICE Macro-Model of memristor 

Basically, there are three different ways available to model the electrical characteristics of 
the memristors. SPICE macro-models, hardware description language (HDL), and C pro- 
gramming. The first, SPICE macro-models, approach is more appropriate since it is more 
readable for most of the readers and available in all SPICE versions. There is also another 
reason for choosing SPICE modeling approach. Regardless of common convergence prob- 
lems in SPICE modeling, we believe it is more appropriate way to describe real device 
operation. Moreover, using the model as a sub-circuit can highly guarantee a reasonably 
high flexibility and scalability features. 



Article submitted to Royal Society 



16 



O. Kavehei et. al. 



A memristor can be realised by connecting an appropriate nonlinear resistor, inductor, 
or capac itor across p ort 2 of an M-R mutator, an M-L mutator, or an M-C mutator, respec- 
tively llChua , 1971 ] . These mutators are nonlinear circuit elements that may be described 
by a SPICE macro-model (i.e. an analog behavioural model of SPICE). The macro circuit 
model realisa tion of a type-1 M-R mutator based on the first realisation of the memristor 
I Chualll97lll is shown in Fig. [8] 




Figure 8. The SPICE macro-model of memristor. Here G, H and S are a Voltage-Controlled Cur- 
rent Source (VCCS), a Current-Controlled Voltage Source (CCVS), and a Switch (V n = -1.9 V 
and Voff = —2 V), respectively. R\ = R2 = I k£l and V c = —2 V. The M-R mutator consists of 
an integrator, a Current-Controlled Voltage Source (CCVS) "H", a differentiator and a Voltage-Con- 
trolled Current Source (VCCS) "G". The nonlinear resistor is also realised with resistors R\, Ri, and 
a switch. Therefore, the branch resistance is 1 kQ. for V < 2 Volt and 2 kQ. for V > 2 Volt. 



In this model the M-R mutator consists of an integrator, a Current-Controlled Voltage 
Source (CCVS) "H", a differentiator and a Voltage-Controlled Current Source (VCCS) 
"G". The nonlinear resistor is also realised with resistors R\, R2, and a switch. Therefore, 
the branch resistance is 1 kQ for V < 2 Volt and 2 kQ for V > 2 Volt. The input voltage 
of port 1, V\, is integrated and connected to port 2 and the nonlinear resistor current, h, is 
sensed with the CCVS "H" and differentiated and converted into current with the VCCS 
"G". 

SPICE simulations with the macro-model of the memristor are shown in Figs.l9land[T0l 
In this particular simulation, a monotonically-increasing and piecewise-linear q-q> function 
is assumed as the memristor characteristic. This function is shown in Fig.|9](b). The simu- 
lated memristor has a value of 1 kQ when the flux is less than 2 Wb, but it becomes 2 kQ 
when the flux is equal or higher than 2 Wb. The critical flux (<p c ) can be varied with the 
turn-on voltage of the switch in the macro-model. Fig. [9] (a) shows the pinched hysteresis 
characteristic of the memristor. The input voltage to the memristor is a ramp with a slope 
of ±1 V/s. When the input voltage ramps up with a slope of ±1 V/s, the memristance 
is 1 kQ and the slope of the current- voltage characteristics is 1 mA/V before the the flux 
reaches to the <p c . But when the flux becomes 2 Wb, the memristance value is changed to 
2 kQ and the slope is now 0.5 mA/V. After the input voltage reaches to the maximum 
point, it ramps down and the slope is maintained, because the memristance is still 2 kQ. 
Fig. [TO] shows the memristor characteristics when a step input voltage is applied. Initially 
the memristance is 1 kQ, so the input current is 1 mA. When the flux reaches to 2 Wb 
(1 V x 2 s), the memristance is 2 kQ and so the input current is now 0.5 mA as predicted. 
The developed macro-model can be used to understand and predict the characteristics of a 
memristor. 

f For further detail about the mutator the reader is referred toChut] Il968ll . 
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Figure 9. The memristor characteristics, (a) The hysteresis characteristics of the memristor. (b) A 
monotonically-increasing and piecewise-linear q-(p function as a basic memristor characteristic, the 
both axes are normalised to their maximum values. The simulated memristor has a value of 1 k£l 
when the flux is less than 2 Wb, but it becomes 2 kfl when the flux is equal or higher than 2 Wb. The 
critical flux (<p c ) can be varied with the turn-on voltage of the switch in the macro-model. The input 
voltage to the memristor is a ramp with a slope of ±1 V/s. 
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Figure 10. The memristor characteristics when a step input voltage is applied, (a) Voltage curve, (b) 
Flux linkage curve, (c) Current curve. At the first point the memristance is 1 kQ., so the input current 
is 1 mA. When the flux reaches to 2 Wb (1 V x 2 s, 2 s = 50 time steps), the memristance is 2 K2 
and so the input current is now 0.5 mA as predicted. 



Now, if a 1 kHz sinusoidal voltage source is connected across the memristor model, 
the flux does not reach to 2 Wb so M = M\ = 1 ki2 and i = 10 mA. Fig. Qjll) shows 
the memristor characteristics when a sinusoidal input voltage is applied. As it is shown in 
Fig. [TlT lI). when the voltage source frequency reduces to 10 Hz, the flux linkage reaches 
to 2 Wb within 30 ms. Based on this result, there are two levels of memristance, M = M\ = 
1 ki2 and then it changes to Mi = 2 ki2. 
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Figure 11. The memristor characteristics when, (I) a 1 kHz sinusoidal voltage is applied. In this case 
the flux does not reach to 2 Wb so M = Mi = 1 kfl and i = 10 mA. (II) When a 10 Hz sinusoidal 
voltage is applied. In this case the flux linkage reaches to 2 Wb within 30 ms, or 15 time steps, (a) 
Voltage curve, (b) Flux linkage curve, (c) Current curve. 

Another interesting study is needed to verify that the model is working properly when 
there is a parallel, series, RM (Resistor-Memristor), LM (Inductor-Memristor), or CM 
(Capacitor-Memristor) network. First of all, let us assume that there are two memristors 
with the same characteristic as shown in Fig. [9] Analysing series and parallel configura- 
tions of these memristors are demonstrated in Figs, \l2\ b) and [T2l a). respectively. In both 
figures, the left I-V curve shows a single memristor. 

The simulation results show that the series and parallel configurations of memristors 
are the same as their resistor counterparts. It means the equivalent memristances, M eq , of 
a two memristor in series and parallel are M eq = 2M and M eq = M/2, respectively, where 
M is memristance of a single memristor. The second step is RM, LM, and CM networks. 
In these cases a 10 V step input voltage has been applied to circuits. As mentioned before 
for a single memristor based on the proposed model, while the flux linkage is less than 
or equal to the critical flux, <f> < <p c , M = Mi = 1 k£2 and when the flux is more than the 
critical flux value, <p > <p £ , M = Mi = 2 k£2. Recall that the critical flux value based on 
the q- (p curve is <p c = 2 Wb. Figs.EtRXEtQ, andQjtL) illustrate RM, CM, and LM 
circuits and their response to the input step voltage, respectively. 

If we assume that the memristance value switches at time T ( /, then for < t < T c (, 
<p < <p £ , and M = M] = 1 kll Therefore, in the RM circuit we have, Vm = V . For 
R = 1 ki2, Vm = 5 V and then im = 5 mA, so 7^ = ^ = 0.4 s. Likewise, when t > 7j, 

<p > <p ( , we have, Vm = V R + M2 = 6.7 V and im = 3.3 mA. Both cases have been verified 
by the simulation results. 

In the LM circuit we have the same situation, so for < t < Tj, q> < q) c , Vm = V = 10 V, 
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Voltage v(t) [V] 

Figure 12. The I-V curves for, (a) two memristors in parallel, (b) two memristor in series, and (c) 
a single memristor. In all of the cases there is no difference between a memristor and equivalent 
resistor in the network. In other words, two memristors in parallel, with the same characteristics, 
form a single memristor with M equ = M q /2 as memristance, and two memristors in series form a 
single memristor with M eqll = 2M q as memristance. 



z'm = 10 mA (R = 1 k£2), and T c [ = 0.2 s. For t > Tj, memristor current is ;'m = 5 mA. 
Memristor's current changing is clearly shown in Fig. [T3lL). The CM circuit simulation 
also verifies a change in time constant from M\C to MjC. 

As another circuit example of using the new memristor model, an op-amp integrator 
has been chosen. The model of an op-amp and circuit configuration is shown in Fig. [14] If 
C = 100 jilC then we have 0.1 s and 0.2 s as the time constant of the circuit at t < Tj and 
t > T c i, respectively. 

It is wor th mentioning that , recen tl y, a few simple SPIC E macro-models have b een 
proposed by iBenderli and Wevl ll2009ll . Biolek et alJ l2009bll and IZhang et all ll2009ll but 
none of them consider the model response with different circuit elements types, which is 
an important step to verify the model correctness. 



4. Interpreting memristor in Electromagnetic Theory 

In his original paper IChual Il97lh presented an argumentation based on electromagnetic 
field theory for the existence of the memristor. His motivation was to interpret the memris- 
tor in terms of the so-called quasi-static expansion of Maxwell's equations. This expansion 
is usually used to give an explanation to the elements of circuit theory within the electro- 
magnetic field theory. 

Chua's argumentation hints that a memristor might exist, it although never proved that 
this device can in fact be realised physically. In the following we describe how Chua argued 
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Figure 13. Step voltage response curves for (R) Resistor-Memristor, (C) Capacitor-Memristor, and 
(L) Inductor-Memristor. (a) Memristor voltage curve, (0- Q 3 ) Flux linkage curve, <p(t). (c) Cur- 
rent curve, i'c(0- The memristance value switches at time 7j, then for < / < 7j, <p < <p c , and 
M = Mi = l k£2. Therefore, in the RM circuit we have, V M = For R = 1 k£l, V M = 5 V and 

then im = 5 mA, so = = 0.4 s. Likewise, when t > 7^, (p > <p c , we have, Vm = V r+m 2 = 6-7 V 
and im = 3.3 mA. In the LM circuit we have the same situation, so for < t < Tj, (p < (p c , 
Vm = V = 10 V, iM = 10 mA (R = 1 k£i), and Tj = 0.2 s. For t > Tj, memristor current is Im = 5 mA. 




20 40 60 80 20 40 60 80 

Time steps Time steps 
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Figure 14. Memristor-op-amp integrator circuit and its response to the input step voltage, (a) Output 
voltage curve, V O ut(0- 0>) F mx linkage curve, <p(t). If C = 100 flC then we have 0.1 s and 0.2 s as 
the time constant of the circuit at t < and t > Tj, respectively. 
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for a memristor from a consideration of quasi-static expansion of Maxwell's equations. To 
consider this expansion, we use Maxwell's equations in their differential form, 



VD = p, (4.1) 
VB = 0, (4.2) 

VxE = — ^, (4.3) 
VxH = J + ^, (4.4) 

where E is electric field intensity (V/m), B is magnetic flux density (Wb/m 2 ), J is elec- 
tric current density (A/m 2 ), p is electric charge density (C/m 3 ), H and D are magnetic 
field intensity (A/m) and electric flux density (C/m 2 ), V- and Vx are divergence and curl 
operators. 

The idea of a quasi-static expansion involves using a process of successive approxima- 
tions for time-varying fields. The process allows us to study electric circuits in which time 
variations of electromagnetic field are slow, which is the case for electric circuits. 

Consider an entire family of electromagnetic fields for which the time rate of change 
is variable. The family of fields can be described by a time-rate parameter a which is time 
rate of change of charge density p. We can express Maxwell's equations in terms of the 
family time, x = at, and the time derivative of B can be written as 

dB dBdx dB 

dt dx dt dx 

Other time derivatives can be expressed similarly. In terms of the family time, Maxwell's 
equations take the form 

VxE = -a^?-, VxH = J+a^, (4.6) 
dx dx 

which allow us to consider different values of the family time x, corresponding to different 
time scales of excitation. Note that in Eqs. 14.61 E. H, D, J, and B are al so functions of 
a and x, along with the position (x,y,z).This allows us i Fano et al. , 1960ll to express, for 
example, E(x,y,z, a, t) as power series in a: 

E(x,y,z,a,t) = E (x,y,z,x) + aEi{x,y,z,x) + a 2 E 2 (x,y,z,x) + (4.7) 
where the zero and first orders are 



E (x,y,z,x) 
Ei(x,y,z,T) 



[E(x,y,z,a,x)] a=0 
dE{x,y,z,a,x) 



da 



a=0 



dx ' 



(4.8) 



Along with these, a similar series of expressions for B, H, J, and D are obtained and can 
be inserted into Eqs. 14. 61 with the assumption that every term in these series is differentiable 
with respect to x, y, z, and x. This assumption permits us to write, for example, 
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V x E = V x E + a(V x Ei) + a 2 (V x E 2 ) + (4.9) 

and, when all terms are collected together on one side, this makes Eqs. l4.6l to take the form 
of a power series in a that is equated to zero. For example, the first equation in Eq. 14.61 
becomes 

V x E + o(V x Ei + + a 2 (V x E 2 + - A + ... = 0, (4.10) 

which must hold for all a. This can be true if the coefficients of all powers of a are 
separately equal to zero. The same applies to the second equation in Eqs. l4.6l and one then 
obtains the so-called n th-order Maxwell's equations, where n = 0, 1,2, ... For instance, the 
zero-order Maxwell's equations are 



VxE = 0, (4.11) 
VxHo = Jo, (4.12) 



and the first-order Maxwell's equations are 



VxEi = (4.13) 
<?D 

VxHi = Ji + V- ( 4 - 14 ) 

OT 

The quasi-static fields are obtained from only the first two terms of the power series 
Eq. 14.101 while ignoring all the remaining terms and by taking a = 1 . In this case we can 
approximate E « E + Ei , D « D + Di , H ss H + Hj , B w B + Bi , and J w J + Ji . 

Circuit theory, along with many other electromagnetic systems, can be explained by the 
zero-order and first-order Maxwell's equations for which one obtains quasi-static fields as 
the solutions. The three classical circuit elements resistor, inductor, and capacitor can then 
be explained as electromagnetic systems whose quasi-static solutions correspond to certain 
combinations of the zero-order and the first-order solutions of Eqs. l4.llTl4.14l 

However, in this quasi-sta tic explanat i on of circuit elements, an interesting possibility 
was unfortunately dismissed IFano et al. . Il960ll as it was thought not to have any corre- 



spondence with an imaginable situation in circuit theory. This is the case when both the 
first-order electric and the first-order magnetic fields are not negligible. Chua argued that it 
is precisely this possibility that provides a hint towards the existence a fourth basic circuit 
device. 

Chua's argumentation goes as follows. Assume there exists a two-terminal device in 
which Di is related to Bi, where these quantities are evaluated instantaneously. If this is 
the case then this device has the following two properties: 

1 . Zero-order electric and magnetic fields are negligible when compared to the first- 
order fields i.e. Eo, Do, Bo, and Jo can be ignored. 

2. The device is made from non-linear material for which the first-order fields become 
related. 
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Assume that the relationships between the first-order fields are expressed as 

Ji = /(Ei), (4.15) 
Bi = #(Hi), (4.16) 
Di = #(Ei), (4.17) 

where and & are one-to-one continuous functions defined over space coordinates 

only. Combining Eq. 14.141 in which we have now Do ~ 0, with Eq. l4.15l gives 

VxH, = ^(Ei). (4.18) 

As the curl operator does not involve time derivatives, and ^ is defined over space coor- 
dinates, Eq. 14. 181 says that the first-order fields Hi and Ei are related. This relation can be 
expressed by assuming a function & and we can write 

Ei=JT(Hi). (4.19) 
Now, Eq. l4.17l can be re-expressed by using Eq. 14. 191 as 

Di = ^o^(Hi), (4.20) 

where o operator is the composition of two (or more) functions. Also, as 53 is a one-to-one 
continuous function, Eq. l4.16l can be re-expressed as 

H^i^Bj). (4.21) 
Inserting from Eq. l4.21l into Eq. l4.20l then gives 

Di =&o^o[S§-\^ l )} =S?(Bi). (4.22) 

Eq. l4.22| predicts that an instantaneous relationship can be established between Di and 
Bi that is realisable in a memristor. This completes Chua's argumentation using Maxwell's 
equations for a quasi-static representation of the electromagnetic field quantities of a mem- 
ristor. 

5. Conclusion 

In this paper, we surveyed key aspects of the memristor as a promising nano-device. 
We also introduced a behavioural and SPICE macro-model for the memristor and re- 
viewed Chua's argumentation for the memristor by performing a quasi-static expansion 
of Maxwell's equations. The SPICE macro-model has been simulated in PSpice and shows 



agreement with the actual memristor presented in lStrukov et al.l 11200811 . The model shows 
expected results in combination with a resistor, capacitor, or inductor. A new op-amp based 
memristor is also presented and tested. 

Nanoelectronics not only deals with the nanometer scale, materials, and devices but 
implies a revolutionary change even in computing algorithms. The Von-Neumann archi- 
tecture is the base architecture of all current computer systems. This architecture will need 
revision for carrying out computation with nano-devices and materials. There are many 
of different components, such as processors, memories, drivers, actuators and so on, but 
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they are poor at mimicking the human brain. Therefore, for the next generation of com- 
puting, choosing a suitable architecture is the first step and requires deep understanding of 
relevant nano-device capabilities. Obviously, different capabilities might create many op- 
portunities as well as challenges. At present, industry has pushed nanoelectronics research 
for highest possible compatibility with current devices and fabrication processes. However, 
the memristor motivates future work in nanoelectronics and nano-computing based on its 
capabilities. 

In this paper we addressed some possible research gaps in the area of the memristor 
and demonstrated that further device and circuit modelling are urgently needed. The cur- 
rent approach to device modelling is to introduce a physical circuit model with a number 
of curve fitting parameters. However, such an approach has the limitation of requiring a 
large number of parameters. Using a non-linear drift model results in more accurate sim- 
ulation at the cost of a much more complicated set of mathematical equations. Initially 
behavioural modelling can be utilised, nonetheless a greater modelling effort is needed to 
accommodate both the defect and process variation issues. An interesting follow up would 
be the development of mapping models based on the memristor to neuronmorphic systems 
that deal with architectural level challenges, such as defect-tolerance and integration into 
current integrated circuit technologies. 
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